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Abstract 


Among the risks facing the plains in Iran is subsidence, which causes many problems in 
agricultural lands, roads, power, and energy transmission lines. In recent years, Sarab plain has 
faced a sharp drop in the level of underground water, which has caused this area to be exposed to 
the risk of subsidence. Therefore, the purpose of this research is to investigate and analyze the 
most important factors involved in creating the risk of subsidence in Sarab plain and to identify 
the susceptible surfaces that are likely to be involved in subsidence in the future, using the multi- 
criteria MARCOS and CODAS algorithms. According to the results of subsidence risk zoning, 
water depth, land use, and slope, respectively, with weight coefficients of 0.194, 0.171, and 0.159, 
are the most important factors involved in creating the risk of subsidence in the studied area. The 
output of the MARCOS method showed that, respectively, 167.50 and 276.09 square kilometers 
of the area of Sarab Plain, and results of applying the CODAS method, showed that 187.13 and 
279.03 square kilometers of the area are in the high-risk and critical category. In addition, the 
map extracted from the MARCOS and CODAS algorithms with the depth of the water level of 
the wells, respectively, have correlation coefficient values of 0.77 and 0.81. A correlation can be 
seen between the output of both methods with the water level map. It seems that the results of this 
study can be of great help to organizational managers and land and soil resource planners for 
protecting and managing water resources and natural hazards and preventing land degradation. 
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جغرافیا و مخاطرات محیطی. سال یازدهم. شمارة چهل و چهارم. زستان ۱۶۰۱ صص ۱2۹-۱۷۱ 
مقاله پژوهشی 


ویژه‌نامه (چالش جهانی فرونشست زمین: مدیریت بحران یا بحران مدیریت) 


پهنه‌بندی $$ خطر فرونشست VE‏ در دشت سراب. با بهره‌گیری از الگوریتم‌های تحلیل چند معیاره 
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حکیده 


یکی از مخاطرات پیش روی دشت‌های کشور. مخاطره فرونشست است که سبب بروز مشکلات و 
معضلات فراوان در زمین‌های کشاورزی. جاده‌هاء خطوط انتقال نیرو و انرژی می‌گردد. دشت سراب نیز در 
طی سال‌های اخیر با افت شدید سطح آب زیرزمینی مواجه بوده که این عامل سبب گردیده تا این منطقه 
در معرض وقوع مخاطره فرونشست قرار گیرد. هدف این پژوهش بررسی و تحلیل مهم‌ترین عوامل دخیل 
در ایجاد خطر فرونشست دشت سراب و مشسخص کردن سطوح مستعد که احتمالاً در آینده درگیر 
فرونشست خواهند شد. با بهره‌گیری از الگوریتم‌های چند معیاره MARCOS‏ و CODAS‏ است. با توجه 


به نتایج حاصل از پهنه‌بندی daz‏ فرونشست؛ معیارهای عمق آب. کاربری اراضی و شیب؛ به ترتیب با 
ضریب وزنی ۰/۱۹۶ ۰/۱۷۱ و ۰/۱۵۹ مهم‌ترین عوامل دخیل در ایجاد خطر فرونشست محدوده مطالعاتی 
می‌باشند. با توجه به خروجی حاصل از روش MARCOS‏ به ترتیب؛ ۰ و ۲۷/۰۹ کیلومترمربع از 
مساحت دشت سراب و طبق نتایج حاصل از به‌کارگیری روش CODAS‏ ۱۸۷/۱۳ و ۲۷۹/۰۳ کیلومترمربع 


Email: s.asghari@uma.ac.i ۰۹۱۰۶۰۲۰۲۵۱ نویسنده مسئول:‎ ۱ 


نحوه ارجاع به این مقاله: 


اصغری سراسکانرود. صیاد؛ پیروزی الناز؛ آقایاری لیلا؛ ۱ مبهنه‌بندی کسترهُ خطر فرونشست زمین در دشت سراب. با بهره گیری 
از الگوریتم‌های تحلیل چند MARCOS. Gs‏ و CODAS‏ جغرافیا و محاطرات محیطی. ۶(۱۱). صص ۱۶۹-۱۷۱ 


https://doi.org/10.22067/geoeh.2023.79331.1291 


o:‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


از مساحت محدوده در dab‏ پرخطر و بحرانی قرار دارد. به‌علاوه. نقشه مستخرج از الگوریتم MARCOS‏ 
و CODAS‏ با عمق سطح آب جاه‌ها به ترتیب؛ دارای مقدار ضریب همبستگی ۷۷/ و ۰/۸۱ می‌باشند و 
بین خروجی هر دو روش با نقشه سطح آب همبستگی دیده می‌شود. به نظر می‌رسد. نتایج حاصل از 
مطالعه حاضر می‌تواند به مدیران سازمانی و برنامه‌ریزان منابع اراضی و خاک در dyes‏ حفاظت و مدیریت 


کلیدواژه‌ها: فرونشست. یهنه‌بندی» abs‏ سراب. الگوریتم‌های تحلیل چند معیاره. 


۱- مقدمه 


فرونشست زمین به‌عنوان پدیه‌ای مورفولوژیک نوعی از تغیبر شکل سطح زمین است که با دگرشکلی عمودی 
و یا حرکت رو به پایین سطح زمین و همچنین نشست تدریجی یا دفعی مواد سطحی همراه است (شریفی کیا و 
همکاران. ۱۳۹۶). فرونشست به fle‏ گوناگونی ازجمله عوامل طبیعی مانند زلزله. آتشفشان, فعالیت‌های گسلی. 
فرونشست ناشی از بالا آمدن سطح دریاء انحلال در سنگ اکسایش» فشردگی رسوبات آلی» توسعه چاله در زمین‌های 


کارستی یا ناشی از فعالیت‌های انسانی» شامل برداشت بی‌رویه سیالات از زمین نظیر آب» نفت و گاز روی می‌دهد 


(شفیعی و همکاران ۱۳۹۹). هرچند محوریت رویداد آن متوجه دو عامل اصلی استخراج آب زیرزمینی و کارستی 
شدن سنگ‌های انحلال‌پذیر است (شریفی LS‏ ۱۳۹۱). فرونشست سبب تخریب کانال‌های آبیاری» شکست خحطوط 
لوله‌های آب و نفت. خحطوط انتقال نیرو و تخریب سازه‌ها و تشدید فرایند زمین‌لغزه و سیلاب‌ها می‌شود و در صورت 
عدم شناخت این پدیده سبب آسیب جدی به زیرساخت‌های اقتصادی می‌گردد (گامبولاتی و همکاران » ۲۰۱۸). 
فرونشست زمین در dl yl‏ سابقه‌ای طولانی دارده در گذشته تنها برحی از استان‌ها ازجمله کرمان و یزد با این مسئله 


مواجه بودند ولی هم اکنون بهره‌برداری بیش‌ازحد از آب‌های زیرزمینی سبب شله تا بسیاری از دشت‌های D sl‏ با 


alton‏ فرونشست مواجه ntl‏ بر آساس آمار اعلام شده در کشور cil gl abl yl‏ سوء ناشی از رخداد فرونشست رقم 

پایینی نیست و به‌سرعت در حال ایجاد و گسترش در دشت‌های مختلف سراسر کشور است که عدم مدیریت و 

کنترل به‌موقع عوامل آن می‌تواند خسارات جانی و مالی جبران‌ناپذیری برجای بگذارد (منتظریون و اصلانی: ۱۳۹۸). 
بايد توجه داشت که از یک‌سی یکی از مهم‌ترین راهکارهای کنترل فرونشسست در مناطق بر اساس تجارب 


جهانی. روش‌های ue‏ سازه‌ای کاهش‌دهنده E‏ تطبیق‌دهنده انیت (آبیدین E‏ اران 10( و از سوی دیگ 


اندازه‌گیری رخداد فرونشست در مقیاس وسیع به دلیل نرخ AS‏ آن چالش‌برانگیز است و نیاز به فناوری‌های پیشرفته 


دارد. آمروزه سیستم اطلاعات جغرافیایی» ابزاری مور برای سیستم‌های پشتیبانی تصمیم گیری تعاملی برای عملیات 


1 Gambolati et al 
2 Abidin et al 


سال یازدهم پهنه بندی گسترهٌ خطر فرونشست زمین در دشت سراب .... Yo‏ 


مدیریت ریسک هستند (بوخیر و همکاران. (TH‏ و همچنین» برای کشف طیف وسیعی از گزینه‌ها از نظر 
درگیری‌های عینی و معیارهای چند lS‏ از روش MCDM‏ استفاده می‌شود (آهر و همکاران ؛ ۲۰۱۳). در این 
روش‌های تصمیم‌گیری برای بهینه‌سازی مدل از معیارهای مختلفی برای افزایش صحت تصمیمات استفاده می‌شود 
(جورجیو و همکاران » ۲۰۱۵؛ عرب عامری و همکاران ۰ ۲۰۱۸). لذا می‌توان گفت. استفاده از سیستم‌های اطلاعات 
جغرافیایی و روش‌های تصمیم‌گیری چندمعیاره با یک رویکرد تلفیقیء می‌تواند باعث تسریع در روند برنامه‌ریزی در 
تشخیص موارد بحرانی و اضطراری شده و منجر به صدور نتایج مناسبی گردد. 


با توجه به اهمیت موضوع در سال‌های اخیرء مطالعات مختلفی در زمینه فرونشست صورت گرفته است. کیانی و 


همکاران (۱۳۹۷) اقدام به بررسی ارتباط فرونشست زمین و افت سطح آب‌های زیرزمینی در شهرستان کرج با 
استفاده از روش تلفیق وزنی در محیط GIS‏ کرده و نتیجه کارشان نشان داده که بین وضعیت توپوگرافی. ضخامت 
سازند و برداشت آب ارتباط تنگاتنگی وجود دارد و بیشترین میزان فرونشست در مناطق با برداشت زیاد آبرفت‌های 
ضخیم و مناطق پست دشت مشاهده می‌شود. منتظریون و اصلانی (۱۳۹۸) با به‌کارگیری سیستم اطلاعات جغرافیایی 
به ارزیابی خطر فرونشست زمین در پهنه استان‌های تهران و البرز پرداختند. یافته‌های پژوهش نشان داد مناطق 
وسیعی در استان تهران ازجمله غرب و جنوب غربی و همچنین دشت ورامین» شهریار کرج» نظرآباد و ساوجبلاغ در 
استان البرز ازجمله مناطق مستعد وقوع فرونشست هستند و آنچه باعث ایجاد خطر فرونشست می‌گردد عدم JES‏ 
و مدیریت مناسب منابع و ذخایر آب‌های زیرزمینی در سطح استان است. ابراهیمی و همکاران (۱۳۹۹)» به پایش 
میزان فرونشست محدوده شهری دشت پاکدشت با روش تداخل‌سنجی پرداختند. نتایج حاصل از مطالعه» در طی بازه 
زمانی یک‌ساله بین ۱۵ تا ۸۵ میلی‌متر فرونشست را نشان می‌دهد و مهم‌ترین و موثرترین عامل در فرونشست را 
هم‌افت سطح آب‌های زیرزمینی عنوان کرده‌اند. فرزین کیا و همکاران OE)‏ پهنه‌بندی خطر فرونشست زمین در 
دشت جوین را با استفاده از مدل تحلیل شبکه‌ای _ فازی موردمطالعه قراردادند. با توجه به نتایج مطالعه حدود 0 
درصد از مساحت محدود و شامل قسمت مرکزی دشت دارای احتمال حطر زیاد تا he‏ زیاد است و تراکم بالای 
چاه‌های عمیق و نیمه عمیق و افت بالای آب‌های زیرزمینی به دلیل برداشست بی‌رویه» به‌عنوان دلیل اصلی استعداد 
بالای فرونشست محدوده مشخص شد. روستایی و همکاران (۱۶۰۱) مناطق مستعد ریسک فرونشست زمین در اثر 
افت سطح آب زیرزمینی با استفاده از روش FUZZY-AHBP‏ در دشت شبستر-صوفیان موردبررسی قرار دادند. طبق 
نتایج مطالعه. مناطق جنوبی و غربی دشت شبسترح صوفیان مستعد فرونشست زیاد تا خیلی زیاد بوده است و این 


1 Aher et al 
2 Georgiou et al 
3 Arab Ameri et al 


oY‏ جغرافیا و مخاطرات محیطی شمارۂ چهارم 


پهنه‌بندی خطر فرونشست زمین در دشت مرکزی استان البرز با بهره‌گیری از تکنیک تداخل‌سنجی راداری و الگوریتم 

تحلیل چند معیاره ARAS‏ پرداختند و نتایج مطالعه OUS‏ داد» فرونشست در محدوده بین صفر تا ۳۰۰ میلی‌متر بوده و 

به ترتیب ۱۳۵/۵۵ و ۱۹۲/۲۸ کیلومترمربع از محدوده دارای احتمال خطر بسیار 5 3b‏ و زیاد ul‏ به‌علاوه مهم‌ترین 

عامل اصلی دخیل در افزایش پتانسیل فرونشست محدوده بهره‌برداری بی‌رویه از آب‌های زیرزمینی و افت سطح آب 

مناطق مسکونی پرداختند و طبق بررسی‌های انجام شده برداشت بی‌رویه منابع زیرزمینی و تجمع سازه‌ها را مهم‌ترین 
بو و "P " , Y " i . ۳ ۳ ba‏ 

علت فرو در مناطق مسکونی بیان کردند. مين و همکاران )14 ۰ در پژوهشی به بررسی فرونشست در 


محدوده شناخته شده ial‏ هوانگ و همکاران؟ (۰۲۰)» به بررسی فرونشست 


مهم‌ترین عامل موثر در فرو 
شهر دژو با روش SBAS‏ پرداختند. بررسی‌ها طی بازه زمانی دو ساله OUS‏ می‌دهد که مرکز فرونشست بخش‌های 
شرقی و شمال‌شرقی با میزان فرونشست £0 میلی‌متری است و این نشست هم رابطه مستقیمی با افت سطح آب‌های 
زیرزمینی دارد. رنجبر و همکاران ‏ (۲۰۲۱) به بررسی حساسیت فرونشست زمین, در شهرستان شهریار با استفاده از 
تکنیک تداخل سنجی و الگوریتم‌های یادگیری ماشین‌بردار پرداختند و به این نتیجه رسیدند در سال‌های اخیر به دلیل 
افزایش جمعیت و برداشت بیشتر از آب‌های زیرزمینی؛ نرخ فرونشس" 
(TTY)‏ در پژوهشی. با استفاده از InSAR SS‏ به مطالعه à jar‏ 235 ت زمین در شهر دهلی هند در طی 
سال‌های ۲۰۱ تا ۲۰۲۰ toh,‏ با توجه به نتایج حاصل از این پژوهش» به‌طور تقریبی» ۰ کیلومترمربع از 


مساحت منطقه در معرض خطر 233 - زمین قرار دارد و گسترش و برداشت بیش‌ازحد آب‌های زیرزمینی» 


افزایش یافته است. گرج و همکاران " 


به‌عنوان علت اصلی فرونشست زمین, در منطقه است. 

دشت سراب. از دشت‌های حاصلخیز استان آذربایجان شرقی است که اقتصاد جوامع ساکن در آن به کشاورزی و 
دامداری وابسته است. در این دشت به دلیل شرایط خشک حاکم بر منطقه و بارش دریافتی کم. کشاورزان علاوه بر 
آب‌های سطحی از آب‌های زیرزمینی به‌عنوان مکمل برای آبیاری مزارع و باغات استفاده می‌کنند (اسفندیاری و 
همکاران؛ ۱۳۹۷). dI‏ در سال‌های اخیر این دشت با مشکل برداشت بی‌رویه آب و درنتیجه فقر سفره‌های آب 
زیرزمینی و فرونشست زمین روبه‌رو شده است. به‌طوری‌که کمیته مخاطرات زمینی سازمان نقشه‌برداری کشور نیز در 


راستای ایفای نقش خود در زمینه GL‏ تغیبرات سطحی پوسته زمین در کشور» طی نامه‌ای نرخ و گستره فرونشست 


1 Gambolati et al 
2 Minh et al 

3 Huang et al 

4 Ranjgar et al 

5 Garg et al 


سال یازدهم پهنه بندی گسترهٌ خطر فرونشست زمین در دشت سراب .... YoY‏ 


این منطقه را به وزارت جهاد کشاورزی, مرکز تحقیقات col‏ مسکن و شهرسازی و فرمانداری شهرستان سراب اعلان 
نموده و نسبت به بروز مخاطرات این پدیده هشدار داده است Cols)‏ سازمان نقشه‌برداری» GL Ert‏ توجه به اهمیت 
موضوع. مطالعاتی نیز در راستای بررسی مسئله فرونشست در محدوده دشت سراب انجام شده است. به‌عنوان‌مثال؛ 
ناموری و صدری‌کیا (۱۳۹۹). با استفاده از تصاویر راداری Sentinel‏ و به کمک تکنیک تداخل‌سنجی راداری نرخ 
فرونشست دشت سراب را ما بین سال‌های ۲۰۱۷ تا ۲۰۲۰ مورد ارزیابی قرار دادند. نتایج نشان داد در دشت سراب 
فرونشست زمین با حداکثر نرخ سالانه حدود ۶ سانتیمتر وجود دارد و فرونشست در بخش‌های جنوبی شهر سراب 
نسبتاً قابل‌توجه است. صدریکیا (VEN)‏ با استفاده از روش تحلیل سری زمانی تداخل‌سنجی راداری در فاصله 
سال‌های ۱۳۸۳ تا ۱۳۹۸ به پایش فرونشست زمین در دشت سراب پرداخته است. طبق نتایج حاصل از این پژوهش 
در دشت سراب بیشینه فرونشستی حدود ۳۰/۵- میلی‌متر در سال در راستای خط دید ماهواره و حدود ELIA‏ 
میلی‌متر در سال در راستای قائم دیده می‌شود و در مناطق دارای فرونشست زیاد. کاهش سطح ایستابی وجود داشته 
و رابطه آماری بین سطح ایستابی و جابه‌جایی‌های سطح زمین معنی دار است و رابطه مستقیم بین دو متغیر وجود 
دارد. همان‌طور که بررسی پیشینه پژوهشی نشان می‌دهد. از یک‌سو مطالعات انجام شده در محدوده دشت سراب؛ 
صرفاً با هدف برآورد نرخ فرونشست بر ah‏ بهره‌گیری از تکنیک تداخل‌سنجی راداری استوار بوده است و تاکنون 
مطالعه‌ای به‌منظور پهنه‌بندی سطح دشت سراب از لحاظ پتانسیل رخداد فرونشست انجام نشده است. از سوی دیگر 
با توجه به اینکه فرونشست. ازجمله مخاطراتی است که دارای طیف وسیعی از گزینه‌ها از نظر درگیری‌های عینی و 
معیارهای چندگانه. هستند. یکی از روش‌های کاربردی و موردتوجه در طی سال‌های اخیر استفاده از روش‌های 
تصمیم‌گیری چند معیاره است؛ لذا هدف این مطالعه شناخت و تحلیل عوامل موثر بر فرونشست دشت سراب و 
مشخص کردن سطوح مستعد فرونشست با ارائه نقشه پهنه‌بندی حطر فرونشست زمین با استفاده از الگوریتم‌های 
تحلیل چند معیاره MARCOS‏ و CODAS‏ و سیستم‌های اطلاعات جغرافیایی است. 


-Y‏ مواد و روش‌ها 

۱-۲- منطقه موردمطالعه 

دشت سراب با مساحت ٩۱۸/4۰‏ کیلومترمربع» بین مدارهای TV‏ درجه و LV‏ دقيقه تا ۳۸ درجه Y‏ دقيقه و ۲۱ 
at‏ شمالی و نصف‌النهارهای ÍV‏ درجه و ۱۰ ARS‏ تا EV‏ درجه و ۵۰ دقیقه شرقی قرار گرفته است (شکل ۱). این 
دشست. در میان توده‌های عظیم آتشفشانی سبلان و رشته‌کوه‌های بزقوش aly‏ شده است و پیوستگی کوهستان سبلان 
و بزقوش» در شرق و شمال» این دشت را به‌صورت چاه بسته درآورده که فقط از گوشه شمال‌غربی باز است. دشت 
سراب شکل کشیده‌ای دارد و کشیدگی آن» در جهت شرقی - غربی می‌باشد. دشت سراب در محور رودخانه 


آجی‌چای و در امتداد شرقی-غربی واقع شده است و سفرة آب‌های آزاد این دشت بیشتر از جریان آبی کوهستان‌های 


Lal > ۱‏ و مخاطرات محیطی شمارة چهارم 


سبلان و بزقوش, تأمین می‌شود. منابع آب سطحی دشت را بیو ک‌چای» تاجیارچای» رازلیق‌چای, وانق‌چای تشکیل 
می‌دهند که از ارتفاعات اطراف سرچشمه می‌گیرند و بعد از روستای اندرآب به‌هم‌پیوسته و با نام آجی‌چای از شرق 
به غرب جریان يافته و درنهایت به دریاچه ارومیه می‌ریزند (تلسچی امیرخیزی و همکاران. ۱۳۹۸). این دشت شش 
درة آبرفتی دارد که رسوبات آن‌ها در سمت پایین‌دست به هم ملحق می‌شوند و دشت سراب را به وجود می‌آورند. 
ازنظر زمین‌ساختی. تشکیل دشت سوبسیدانس سراب وابسته به فاز تکتونیکی میوپلیوسن است و این حوضه 
اغلب از رسوب‌های میوسن. پلیوسن کواترنر با ضخامت زياد پوشیده شده است. گسل‌های متعددی در 
شکل‌گیری آن مژثر بوده‌اند که در این میان گسل‌های دامنه شمالی بزقوش و دامنه جنوب سبلان, نقش اساسی در 
شکل‌گیری آن داشته‌اند. سازندهای زمین‌شناسی که در منطقه برونزد دارند. به ترتیب سنی از قدیم به جدید شامل؛ 
کنگلومرا و ماسه‌سنگ به سن پلیومیوسن NBC‏ مارن ژیپسی ماسه‌سنگ دار به سن میوسن Mmg‏ ولکانیک‌های 
پلیوسن ۲1۷ کنگلومرا و ولکانیک‌های پلیوکواترنر ۳0۷ و PLAC‏ و پادگانه‌های آبرفتی Qt‏ به سن کواترنر است. نوع 


اقلیم دشت نیز بر اساس روش آمبرژه نیمه‌عشک سرد است. 
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شکل ۱- نقشه موقعیت محدوده موردمطالعه 
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Y- Y‏ - داده‌ها 3 ابزارهای مورد استفاده 


پژوهش حاضر از نوع کاربردی بوده و روش تحقیق آن تحلیلی مبتنی بر تلفیق آنالیز داده‌ها؛ سیستم اطلاعات 
جغرافیایی و استفاده از فنون تحلیل چند معیاره اسست. جهت پردازش تصاویر و تجزیه‌وتحلیل داده‌ها نیز از 
نرم‌افزارهای Idrisi ENVI, Ecognition, Arc GIS‏ و Excel‏ بهره‌گیری شده است. جهت پهنه‌بندی, حطر 
tij‏ ت» در مرحله جمع‌آوری داده‌هاء ابتدا باید مؤلفه‌های تأثیرگذار مشخص شده و بر مبنای آن‌ها داده‌ها و 
اطلاعات موردنیاز گردآوری و طبقه‌بندی شوند. بدین منظور, در این مطالعه ابتدا عوامل موثر (شامل: شیب لیتولوژی؛ 
کاربری اراضی, بارش فاصله از شهر و روستاء فاصله از رودخانه. فاصله از گسل و افت سطح آب زیرزمینی) با 
بررسی نمونه مطالعات انجام شده در رابطه با خطر فرونشست که توسط دیگر پژوهشگران صورت پذیرفته و عوامل 
مؤثر در ایجاد فرونشست را بیان داشته‌اند ((مانندء کیانی و همکاران (۱۳۹۷» منتظریون و اصلانی YA)‏ ابراهیمی و 
همکاران YA)‏ فرزین LS‏ و همکاران (VEO)‏ اصغری سراسکانرود و همکاران (VEN)‏ روستایی و همکاران 


ER)‏ گامبولاتی و همکاران (PHVA)‏ تانگ مین و همکاران (۲۰۱۹)» هوانگ و همکاران Y)‏ رنجبر و همکاران 
(۲۰۲۱) و گرج و همکاران (۲۰۲۲)) و با توجه به شرایط طبیعی و انسانی منطقه» شسناسایی گردید. در مرحله بعد 
لایه‌های اطلاعاتی مربوط به هر یک از عوامل» در محیط سیستم اطلاعات جغرافیایی تهیه گردید. لایه‌های اطلاعاتی 
شبکه wal yf‏ با استفاده از نقشة رودخانه‌های استان آذربایجان شرقی, تهیه گردید. لایه‌های شیب نیز با استفاده از مدل 
رقومی ارتفاعی ۱۳/۵ متر» دانلود شده از سایت (vertex.daac.asf.alaska.edu)‏ تهیه شد. لایه‌های اطلاعاتی مربوط 
به لیتولوژی (مقاومت سنگ‌ها) و گسل‌هاء با رقومی‌سازی از روی نقشه زمین‌شناسی استان؛ با مقیاس ۱:۱۰۰۰۰۰ تهیه 
شد. برای استخراج کاربری اراضی محدوده موردمطالعه, از تصاویر لندست هشت» سنجنده COLT)‏ اخذ شده از 
سایت cearthexplorer.usgs.gov)‏ استفاده شده است. در این راستاء ابتدا تصحیحات هندسی و اتمسفری بر روی 
تصاویر با استفاده از روش Flaash‏ در نرم‌افزار Envi‏ انجام شد. سپس طبقه‌بندی با روش شیءگرا و الگوریتم 
نزدیک‌ترین همسایگی در نرم‌افزار Ecognition Developer‏ صورت گرفت و نتایج حاصله از طبقه‌بندی کاربری‌ها؛ 
چه ازنظر تک‌تک کاربری‌ها و چه ازنظر مجموع صحت و آمار US‏ از صحت قابل قبولی (بزرگ‌تر از ۸۵ درصد) در 
ارتباط با اطلاعات تولیدشده برخوردار می‌باشد. به‌منظور پایش سطح آب‌های زیرزمینی منطقه موردمطالعه. از اطلاعات 
مربوط به ۳6 چاه مشاهده‌ای سطح تراز پیزومتری دشت سراب Jer)‏ شده از سازمان آب منطقه‌ای استان آذربایجان 
شرقی) استفاده گردید و سپس با توجه به اینکه تخمینگر کریجینگ بهترین تخمینگر خطی cab‏ است (قربانی و 
همکاران. AYAV‏ و با نظر به مقادیر ضریب تعیین UL‏ میزان خطا و واریانس تخمینی کمتر این روش در تهیه نقشه 
سطح آب نسبت به دیگر روش‌ها که (کماسی و همکاران؛ Y‏ (اسفندیاری و همکاران, (WV‏ (اصغری 
سراسکانرود و همکاران» VEY‏ نیز در مطالعات خود بدان اشاره داشته‌اند» درنهایت با استفاده از روش کریجینگ» 
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نقشه عمق آب دشت سراب استخراج گردید. سپس با اجرای مراحل تکنیک مارکوس و کوداس نسبت به تهیه نقشه 
پهنه‌بندی خطر فرونشست اقدام شد. 

۳-۲- مراحل angi‏ نقشه پهنه‌بندی خطر فرونشست با استفاده از روش مارکوس CMARCOS)‏ 

۱- تعبین معیارها و گزینه‌های پژوهش: اولین گام در این روش تعیین عوامل و گزینه‌های پژوهش اسست. -Y‏ 
تشکیل ماتریس تصمیم: ماتریس تصمیم در این روش نیز به‌صورت معیار-گزینه است؛ یعنی یک ماتریسی که 
ستون‌های ol‏ را معیارهای مساله و سطرها را گزینه‌ها تشکیل می‌دهند و هر سلول نیز در واقع امتیاز هر گزینه نسبت 
به هر معیار است. ۳- تعیین ایده آل و ضد ایده آل: در این بخش بر اساس روابط Y‏ ۲ مقادیر sul‏ آل (AI)‏ و ضد 
ایده آل (AAT)‏ مشخص می‌شود. در روابط زیر عبارت 8 به معنی معیارهایی که جنبه سود و عبارت € به معنی 


معیارهای که جنبه هزینه دارند. 


AI = max xj if j € Bandminx;; if j € C Y رابطه‎ 


AAI = min xj; if j € Band maxx; if j € C Y رابطه‎ 


۶- نرمال کردن ماتریس تصمیم یا استانداردسازی: در پژوهش حاضر» جهت دستیابی به الگوی پهنه‌بندی مکانی 
خطر زمینلغزش. مرحله استانداردسازی با توجه به تابع عضویت فازی صورت پذیرفت. ۵- وزن‌دار کردن ماتریس 
نرمال: در این پژوهش. جهت وزن‌دهی عوامل از روش CRITIC‏ استفاده شده است. =I‏ محاسبه مطلوبیت هر گزینه: 
در این مرحله» بر اساس روابط ۳و 4 درجه مطلوبیت ایده K+ JE‏ و ضد ایده آل K-‏ گزینه‌ها محاسبه می‌شود. 


Kj = 2 Y رابطه‎ 
Si 
BE = Saai : sab 


del; 22‏ فوق. PE ET Si=(i=1,2,3,..m)‏ هر pum‏ در ماتریس وزن‌دار می‌باشد که از رابطه <ð‏ به دست 


رابطه ۵ Xja Vij‏ = رک 


۷-تعیین عملکرد نهایی و رتبه‌بندی گزینه‌ها: در این بخش با استفاده از رابطه A‏ عملکرد مطلوب هر گزینه 


محاسبه می‌شود. 


1 Measurement Alternatives and Ranking according to Compromise Solution 
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kf +k; 
PE) ۳ rr), A) 


1+ 


٩ رابطه‎ 


در رابطه فوق‌الذکر Shoe d(k-)‏ > مطلوبیت ضد ایده- آل و (+10 عملکرد مطلوب ایده آل برای هر گزینه است 
که از روابط ۷ و A‏ محاسبه می‌شود. سپس بر اساس اعداد به دست آمده از هر گزینه رتبه‌بندی انجام می‌شود. 


هرکدام از گزینه‌ها مقدار 1000 بزرگ‌تری داشت رتبه بهتری دارد (جهانگیری۱۳۹۹؛ جهانگیری» ۱4۰۰). 


2 kt 
f (kj ) = kfekp 


۳ k; 


V رابطه‎ 


۲-- مراحل تهیه نقشه پهنه‌بندی خطر فرونشست با استفاده از روش کوداس CCODAS)‏ 

۱- تعبین معیارها و گزینه‌های پژوهش: اولین گام در این روش تعبین عوامل و گزینه‌های پژوهش اسست. ۲- 
دومین گام در این روش تشکیل ماتریس تصمیم است. در ماتریس تصمیم روش کوداس معیارهای OF‏ ستون‌ها و 
استاندارسازی و ارزش گذاری به‌صورت als‏ با استفاده از روش فازی انجام شده -í Se cul‏ تشکیل ماتریس نرمال 
وزن‌دار: باید وزن معیارها را در ماتریس نرمال ضرب کرد. در این پژوهش جهت وزن‌دار کردن ماترس نرمال» از 
روش CRITIC‏ استفاده شده است. ۵- محاسبه فواصل اقلیدسی (Euclidean distance)‏ و تاکسی (Taxicab‏ 
ot! jidistances)‏ آل منفی: این فواصل از روابط ٩و‏ ۱۰ به دست می‌آیند. در این روابط nj‏ ایده آل منفی معیارها 


seal 


E = |XjaGi-ns)? ٩ رابطه‎ 
T; = رای(‎ - ns] ۱۰ رابطه‎ 


71 ایجاد ماتریس ارزیابی نسبی: در این گام باید ماتریس ارزیاپی نسبی را با استفاده از رابطه ۱ ایجاد شود که در 
این رابطه ۴ نشان‌دهنده یک تابع آستانه برای تشخیص برابری فاصله اقلیدسی دو گزینه است. 


hig = (Ej — Ex) + (CE; — Ex) x (Tj - Ty) ۱۱ رابطه‎ 


1 Combinative Distance-based Assessment 


۱5۸ جغرافیا و مخاطرات محیطی b‏ چهارم 


۷- در این گام با جمع مقادیر hic‏ گزینه‌هاء می‌توان آن‌ها را رتبه‌بندی نمود هر چه مقدار Hi‏ بزرگ‌تر باشد گزینه 
رتبه بهتری دارد (متیو و ساهو (PR‏ 

-0-Y‏ وزن‌دهی "CRITIC‏ (اهمیت معیارهاء از طریق همبستگی بین معیارها) 

در این روش. داده‌ها بر اساس ميزان تداخل و تضاد موجود بین عوامل يا معیارها مورد تجزیه‌وتحلیل قرار 
فی AS‏ در این روش پس از محاسبهة انحراف معیار معیارهای موردبررسی. ماتریس متفارنی به ابعاد MXM‏ ایجاد 
می گردد که شامل ضرایب همبستگی بین بردارهای تشکیل شده است. با تعیین پارامترهای WL‏ تضاد موجود بین معیار 
Lj‏ معیارهای دیگر از روی Y Abels‏ محاسبه می‌شود. 


Ck-YR - rj) ۱۲ رابطه‎ 


که در آن Cjk‏ معرف مجموع تضاد معیار Lj‏ معیارهای »1 است که از 1 = 6 شروع شده و تا m‏ 2 1 ادامه دارد و 
)زه همبستگی بین دو معیار K‏ و زرا نشان می‌دهد. میزان اطلاعات عامل زرا می‌توان با استفاده از رابطة AY‏ محاسبه 
نمود. 


Cj= زو‎ ELA - rjl) ۱۳ رابطه‎ 


که در آن Cj‏ معرف میزان اطْلاعات معیار زو انحراف معیار در polis‏ مربوط به عامل یا معیار زرا نان می‌دهد. 
با توجه به روابط یادشده معیارهایی که دارای C]‏ بیشتری باشند وزن زیادی به خود اختصاص خواهند داد. وزن هر 
عامل مانند ژاز رابطة ۱۶ تعیین می‌گردد. 


Cj 
wry Ck 


Wj = ۱۶ رابطه‎ 


که در آن Wj‏ معرف وزن معیار زو Ck‏ نشانگر میزان اطلاعات مجموع معیارهای K‏ است که از 1 = ۸ شروع شده 
تمامی معیارها به دست می‌آید انیت (علی‌نداد و خحلیلی ‏ ۹ در راستای انجام مراحل وزن‌دهی CSS S‏ ابتدا در 
محیط ادریسی انحراف معیار و همبستگی میان لایه‌های اطلاعاتی مربوطه به هر یک از عوامل, با استفاده از تابع 
5 به دست آمد و در مرحله بعد با توجه به روابط مربوطه نسبت به محاسبه ماتریس تضاد. Ol xa‏ اطلاعات 
و وزن نهایی در اکسل» اقدام گردید. 


1 Mathew & Sahu 
2 CRiteria Importance Through Intercriteria Correlation 
3 Alinezhad & Khalili 
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—V- Y‏ استانداردسازی فازی 
در پیاده‌سازی تمامی الگوریتم‌های مورد استفاده در پژوهش حاضر و جهت دستیابی به الگوی پهنه‌بندی مکانی 
خطر فرونشست در دشت سراب مرحله استانداردسازی, با توجه به تابع عضویت فازی انجام شده است (جدول (). 
استاندارد نمودن داده‌ها به معنی همسان‌کردن دامنۀ تغییرات داده‌ها بین صفر و یک و یا یک als‏ مشخص دیگر 
است. در مجموعه‌های فازی» بیشترین ارزش یعنی مقدار یک به حداکثر عضویت و کمترین ارزش یعنی صفر به 


حداقل عضویت در مجموعه تعلّق می‌گیرد (حسینی | (M‏ 


جدول ۱- نوع توابع فازی‌ساز لایه‌های موضوعی تأثیر گذار بر خطر فرونشست محدوده موردمطالعه 


متغیر توضیحات نوع ab‏ فازی 

شیب با کاهش مقدار شیب. خطر فرونشست افزایش می‌یابد. تابع خطی کاهشی 

کاربری اراضی اختصاص کدها بر اساس اهمیت آن در وقوع خطر فرونشست و سپس فازی کردن. تابع بزرگ فازی 
فاصله از شهر و روستا در مناطق مجاور از شهر و روستاه پتانسیل وقوع مخاطره فرونشست افزایش می‌یابد. تابع نزدیک فازی 
فاصله از گسل در مناطق مجاور خطوط گسل, پتانسیل وقوع مخاطره فرونشست افزايش می‌یابد. تابع نزدیک فازی 
لیتولوژی اختصاص کدها یر اساس اهمیت آن در وقوع خطر فرونشست و سپس فازی کردن. تابع بزرگ فازی 
فاصله از رودخانه در مناطق مجاور رودخانه پتانسیل وقوع مخاطره فرونشست کاهش می‌یابد. be al‏ کاهشی 
افت سطح آب با افزایش مقدار افت سطح آب» خطر فرونشست افزایش می‌یابد. تابع حطی افزایشی 


۳- نتایج و بحث 
(جدول ۲) و با اعمال مراحل مختلف تکنیک‌های مارکوس و کوداس, نقشه‌های پهنه‌بندی خطر فرونشست در ۵ 
Aib‏ بسیار پرخطر تا بسیار کم خطر به دست آمد (شکل ۲. با am‏ به نتایج حاصله» به ترتیب معیارهای افت سطح 


آب. کاربری اراضیء شیب و زمین‌شناسی بیش‌ترین وزن و اهمیت را دریافت کردند و با توجه به اینکه در وزن‌دهی 
کرتیک. وزن معیارهاء با توجه به نقشه‌های معیار و میزان همبستگی» تضاد و انحراف sles‏ بین لایه‌های اطلاعاتی 
مورداستفاده» تعیین می‌گردد. استفاده از این روش؛ در وزن‌دهی معیارها در پژوهش حاضر می‌تواند معضل استقلال 
صفات از یکدیگر که به هنگام مقایسۂ زوجی در چارچوب روش فرایند تحلیل سلسله مراتبی و فرآیند تحلیل 


شبکه‌ای وجود دارد را برطرف نماید. 
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جغرافیا و مخاطرات محیطی 
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lege y‏ آتنفتالی آندز بنی داسینی | آندزبت و بازالت ولگانیگی بر 
ca jail‏ نیمه ولکانیکی و داسبت EEB‏ آندزیت و توف بای 
heads‏ نفو دی or‏ توف آنهزینی GE‏ 


کتکلومرای درشت AD‏ با la‏ لابه مواد آتشغشانی آواری موضعی | 
رسوبات تراس های آبرفتی و جوان رودخانه ای HEH‏ 
عارن قرمز و قهوه ای هارن ژییی ۵ ار با مبان لابه ماسه سنگ | 
رسوبات تراس gle‏ آبرفتی و رودخانه ای ent‏ 153 
ماسه سنگ و کتکلوتر! ۰ ۱ 
هارن فرمز- عارن ز یپس دار — کتکلومرا و عاسه سنگ |" 


شکل ۲- نقشه معیارهای مطرح در پهنه‌بندی خطر فرونشست دشت سراب 


جدول ۲- مجموع تضاد. انحراف معیان میزان اطلاعات و وزن نهایی معیارهای مطرح 


معیار مجموع تضاد ate cil sl‏ میزان اطلاعات وزن نهایی 
زمین‌شناسی QM NY “YoY WANE‏ 
کاربری اراضی WAVE‏ 00“ ۱/۳۳ ۷۱ 
فاصله از گسل MS UN ۳۵ vvv‏ 
شیب EG: ۳۳ NEN‏ ۰/۹ 
فاصله از رودخانه ۳/۸۷۸ ۳۹۵ WA‏ 7 
سطح آب ‘ov Y/YYA‏ ۱1/10۹ ۰/۱۹۶ 
فاصله از شهر و روستا M NM ۰/۳۹۸ WAVY‏ 


با توجه به حروجی حاصل از به کارگیری روش مارکوس:۱۹۷/۵۰ کیلومترمربع از مساحت محدوده در کلاس 
بسیار پرخطر و ۲۷۹/۰۹ کیلومترمربع در کلاس پرخطر قرار دارد و طبق نتایج حاصل از به‌کارگیری روش کوداس؛ 
۳ و ۲۷۹/۰۳ کیلومترمربع از مساحت محدوده در AE‏ بسیار پرخطر و پرخطر قرار دارد (جدول AY.‏ به‌طور 
کلی؛ با توجه به نتایج هر دو روش مورداستفاده» بخش اعظمی از کلاس بسیار پرخطر و پرخطر از لحاظ احتمال 


وقوع فرونشست با راستای شرقیسفربی و در قسمت مرکزی دشت. قرار گرفته است. 
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شکل ۳- نقشه پهنه‌بندی خطر فرونشست دشت سراب الف: با استفاده از الگوریتم MARCOS‏ ب: با استفاده از 
الگوریتم CODAS‏ 


جدول -Y‏ اطلاعات طبقات خطر فرونشست دشت سراب با استفاده از الگوریتم‌های تصمیم pS‏ 5( چند معیاره 
MARCOS ; CODAS‏ 


طبقة خطر بسیار پرخطر پرخطر خطر متوسط کم خطر بسیار کم خطر 
مساحت به کیلومترمریع ۱/۰ T ۱5۳/۳۸ ۳۹/۸ wed‏ 
MARCOS‏ 
VV ۱3/۰ wir ۳/۰۹ VE die 3d a laa‏ 
مساحت به کیلومترمریع ۱۸/۱۳ wary‏ ۳۳۹/۰۹ ۱۳/۲۸ ۷۹/۹۰ 
CODAS‏ 
سات NV: ۱0/٦ Y£AE ۳/۳۸ ۳/۳۸ die ibis‏ 


می‌باشد و نوع کاربری زمین» خصوصیات سطحی زمین را تحت تأثیر خود قرارداده و موجب تغییر رفتار مواد آن در 
بسیار پرخطر و پرخطر معرفی شده توسط هر دو الگوریتم مارکوس و کوداس؛ OUS‏ می‌دهد؛ از لحاظ معیار کاربری 
را به خود اختصاص داده است. دشت سراب شهر سراب را احاطه کرده و دارای کشتزارها و مراتع وسیع است و 
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شیبی است که شیب آن به کمتر از ۵ درصد می‌رسد)» وجود رسوبات آبرفتی مناسب در کنار عوامل اقتصادی منطقه 
سبب شده است تا فعالیت‌های کشاورزی این منطقه با گسترش زیادی همراه باشد (شکل ۶) به‌طوری‌که بر اساس 
نقشه کاربری اراضی محدوده مطالعاتی» به ترتیب؛ ۲۷۵/۸۶ و 44۲/۲۲ کیلومترمربع از مساحت دشت به اراضی 
زراعت آبی و دیم. احتصاص داده شده است و سراب را به یکی از قطب‌های کشاورزی و دامداری استان آذربایجان 
شرقی تبدیل کرده است و حتی زمین‌هایی که قبلاً در آن‌ها کشت صورت نگرفته بود با ورود ماشینآلات جدید 
کشاورزی به منطقه شخم زده شده است که همین عامل سبب بهره‌برداری بیش ازحد از سفره آب‌های زیرزمینی و 
افزایش پتانسیل رخداد خطر فرونشست. در طی سال‌های اخیر شده است. 


شکل -E‏ تصاویری از پراکنش اراضی کشاورزی واقع در دشت سراب 


می‌توان بیان داشت که عوامل انسان‌ساخت مانند؛ زمین‌های کشاورزی, خطوط انتقال نیرو و عملیات حفاری» 
سکونت گاه‌های شهری و روستایی و شبکه‌های حمل‌ونقل جاده‌ای با افزایش مصرف آب‌های سطحی, زیرزمینی؛ 
تحمیل بار وارده بر زمین در اثر احداث ساختمان‌های بزرگ و تراکم زیرساخت‌ها میزان آسیب‌پذیری این واحد را 
افزایش داده و سبب تشدید فرونشست می‌شوند و در دشت سراب بخش اعظمی از آب موردنیاز جهت شرب نقاط 
شهری و روستایی, از منابع آب زیرزمینی تأمین می‌گردد و این امر نیز به نوبه خود در کاهش سطح آب به دنبال آن, 
در افزایش پتانسیل dam‏ فرونشسست در محدوده موردمطالعه دخیل می‌باشد. به‌طور dS‏ فاصله از مناطق شهری و 
روستایی در محدوده مطالعاتی بین صفر تا ۱۰۰۰ متر می‌باشد. از لحاظ معیار فاصله از مناطق شهری و روستایی بر 
اساس خروجی حاصل از روش مارکوس, مناطق بسیار پرخطر در فاصله صفر تا ۲۵۰۰ متری و Gb‏ نتایج حاصل از 


به‌کارگیری روش کوداس» در فاصله صفر تا ۳۰۰۰ متری از نقاط شهری و روستایی قرار دارند و در این مناطق به 


Lal > ۳4‏ و مخاطرات محیطی شمارة چهارم 


سبب تأمین نیازهای آبی جامعه سکونتگاهی و مناطق اطراف و همچنین تأمين نیازهای آبی بخش کش.ورزی» 
دست‌اندازی به منابع آبی بیشتر بوده و استحصال آب‌های زیرزمینی بسیار بالاست به همین دلیل ميزان افت نیز UL‏ 
بوده است. نمونه‌هایی از بهره‌برداری‌های بی‌رویه از منابع آب زیرزمینی به‌صورت لوله‌کشی و پمپاژ آب چاه‌های حفر 


شده و تأمین نبازهای آبی اراضی کشاورزی واقع در QS‏ سراب؛ در شکل <Ô‏ نمایش داده شده است. 


شکل ۵- نمونه‌هایی از استحصال آب‌های زیرزمینی و تأمین نیازهای آبی اراضی کشاورزی موجود در دشت 


سراب 


اعمال فشار بیش ازحد بر آبخوان ممکن است موجب بروز یدید فرونشست شود؛ ازاین‌رو به بررسی وضعیت 
سفرةٌ آب زیرزمینی در محدودۀ موردمطالعه پرداخته شد. با توجه به نقشه افت سطح آب مستخرح شده از داده‌های 
چاه‌های پیزومتری دشت سراب بیشترین افت سطح آب در بخش میانی دشت (از عمق ۱۲ متر به ۰/۳۰ متر) مشاهده 
می‌شود که احتمال وقوع فروشست G3‏ در این قسمت از محدوده افزایش می‌یابد. ارتباط نقشه‌های نهایی حاصل از 
پهنه‌بندی خطر فرونشست با نقشه معیار سطح آب محدوده موردبررسی قرار گرفت و طبق نتایج حاصله نقشه 
مستخرج از الگوریتم MARCOS‏ با عمق سطح آب چاه‌ها دارای مقدار ضریب همبستگی ۰۸۷۷ می‌باشد (شکل A‏ 
و نقشه پهنه‌بندی فرونشست به دست آمده از الگوریتم CODAS‏ با نقشه عمق سطح آب چاه‌ها؛ ضریب همبستگی 
۷ را نشان می‌دهد (شکل ۷. می‌توان بیان داشت که بین خروجی هر دو روش با نقشه سطح آب همبستگی دیده 
می‌شود. صدری کیا COE Y)‏ نیز در پژوهش خود؛ به این نتیجه دست يافته است که در دشت سراب. در مناطق دارای 
فرونشست ob‏ کاهش چندم‌تری سطح ایستابی وجود داشته و بین سطح ایستابی و جابه‌جایی‌های سطح زمین رابطه 
آماری معنی‌دار وجود دارد و رابطه بین دو متغیر مستقیم است. به‌علاوه ضریب تعیین روش WAY CODAS‏ درصد 


و ضریب تعیین روش MARCOS‏ ۸ درصد می‌باشد؛ به عبارت Gb So‏ نتیجه حاصل از روش کوداس 


سال یازدهم پهنه بندی کته خحطر فرونشست زمین در دشت سراب مب ۱1۵ 


نقشه سطح آب. مربوط به معیار سطح آب است و ۳۳/۱۹ درصد از علت رخداد فرونشست دشت سراب به دیگر 
عوامل مربوط می‌باشد. همچنین؛ Gb‏ مقایسه انجام شده بین نقشه نهایی پهنه‌بندی حاصل از روش مارکوس با نقشه 
سطح آب. معیار سطح آب. ۵۹/۵۸ درصد از علت وقوع فرونشست را OLS‏ می‌دهد و 2۰/۶۲ درصد از احتمال 
فرونشست محدوده با دیگر معیارها در ارتباط است. با نظر به نقش و اهمیت معیار افت سطح آب و با توجه به اینکه 


معیار سطح آب بیش‌ترین وزن و اهمیت را در پتانسیل وقوع فرونشست در سطح دشت سراب دارد. می‌توان بیان 
کرد که میزان موفقیت روش CODAS‏ با توجه به مقایسه خروجی نهایی مدل با معیار سطح آب. از 


روش MARCOS‏ بیشتر است: 


۱ Regress Cal of | 


^ Ye 0040231 +" 0,862532 = t= 0771912 رت‎ EES 
X axia: Water level drop 
Y axis: Subsidence soning - MARCOS 


Comff. of Dat. - 59.580 * 
Std. Dev. of X - 0.357583 
Bud. Dev. of Y - 0.317992 
6.8. of Estimate - 0.202167 
Std. Error of Beta - 0.000976 
t Stat for r or Reta = 884 , 
v 8vav for Bara <> ) = 80.800712 
Sample Size in) = 6 
Apparent df = 530134 


f Show Regression Line 


Save to Clipboard | Close 
Print | Help | 


low frequency high frequency 


شکل -I‏ رابطه همبستگی بین نقشه عمق سطح آب و خروجی حاصل از الگوریتم MARCOS‏ 


Regression Darameters: 


X axis: Water level drop 
Y axis: Subsidence zoning - CODAS 


Coeff. of Dat. - 66.81 4 
Std. Dev. of X E 0.857583 
Std. Dev. of Y - 0.374483 
S.E. of Estimate - 0.215746 
Std, Error of Reta - 0.001041 
© Stat for r or Beca = 1032.998329 
t Stat for leta ir 1 = 72.503277 
Sample fiza (n) - 6 
Apparent df = 630134 


IF Show Hegreseion Line 


Save to Clipboard | Close | 
Print | Help 


low frequen high Irequenc: 


شکل ۷- رابطه همبستگی بین نقشه عمق سطح آب و خروجی حاصل از الگوریتم CODAS‏ 


"n‏ جغرافیا و مخاطرات محبطی شمارة چهارم 


ویژگی‌های سنگ‌شناسی به‌طور گسترده بر ویژگی‌های فیزیکی» مانند؛ استحکام و نفوذپذیری مواد سطحی و 
زیرسطحی و احتمال خطر فرونشست تأثیر می‌گذارد. از منظر سازندهای زمین‌شناسی» به‌طور کلی شهر سراب و 
روستاهای موجود در محدوده دشت سراب. بر روی رسوبات آبرفتی شکل گرفته‌اند و به‌طور عمده La‏ مساحت 
دشت سراب از واحدهای کواترنری (رسوبات پادگانه‌هاای آبرفتی و رودخانه‌ای قدیمی به مساحت We‏ 
کیلومترمربع رسوبات پادگانه‌های آبرفتی و رودخانه‌ای جدید به مساحت ۵۷/۰۸ کیلومتر) و آندزیت و بازالت 
آتشفشانی به مساحت ۶۷/۵۰ کیلومترمربع و مارن قرمز و قهوه‌ای مارن ژیپس‌دار با میان لایه‌های ماسه سنگ به 
مساحت ۳۳/۶۷ کیلومترمربع. است. بررسی مناطق بسیار پرخطر و پرخطر معرفی شده توسط هر دو الگوریتم نیز 
نشان می‌دهد؛ از LES‏ معیار زمین‌شناسی, نقاط دارای احتمال حطر در رسوبات پادگان هه ای آبرفتی و رودخانه‌ای 
قدیمی و مارن قرمز و قهوه‌ای» مارن ژیپس‌دار با ke‏ لایه‌های ماسه سنگ قراردارند و در سایر بخش‌ها به واسطه 
شیب نسبی بیش‌تر و نبود رسوبات آبرفتی زمینه برای وقوع فرونشست بسیار کمتر است. مقدار شیب در محدوده 
موردمطالعه. بین صفر تا YO‏ درصد می‌باشد. محدوده شیب ۰-۲ کمترین میزان شیب در محدوده منطقه می‌باشد که 
این شیب نقش موثری در جمع‌آوری منابع آب زیرزمینی فراهم می‌آورد و هر چه میزان شیب افزایش می‌یابد از میزان 
تمرکز چاه‌ها کاسته و میزان فرونشست زمین هم کاهش می‌یابد. بر اساس نقشه پهنه‌بندی حاصل از هر دو الگوریتم 
مورد استفاده در این پژوهش. از لحاظ معیار شیب مناطق دارای احتمال خطر زیاد و بسیار زیاده بین شیب‌های صفر 
تا ۵ درصد قرار دارند و به‌طور عمده بیشترین مقدار پیکسل‌های با لحتمال خط دارای شیب کمتر از Y‏ درصد 
هستند. 

به دلیل اثر گذاشتن فعالیت گسل‌ها (نظیر زلزله) بر روی فرونشست زمین» بررسی گسل‌های منطقه ضرورت 
می‌یابد. گسل‌ها به‌نوبه خود نقش مؤثر و مهمی در میزان توسعه فرونشست در سطح دشت دارد. به‌گونه‌ای که هر چه 
میزان فاصله از حطوط گسل بیش‌تر باشد نشان‌دهنده این است که منطقه تناسب کمتری برای duel‏ فرونشست 
زمین را دارا بوده و هر چه میزان فاصله به خطوط گسل نزدیک‌تر شود احتمال وقوع فرونشست بیش‌تر می‌شود. 
گسل‌های فراوانی در جهت مختلف در واحدهای سنگی اطراف دشت سراب دیده می‌شود. ولی مهم‌ترین 
آن‌هاء گسل میانی دشت سراب با جهت تقریبی شرقی-غربی است و تقریباً به موازات رودخانه (gle el‏ امنداد 
یافته است. این گسل از نوع فعال است و نشانه فعال بودن OF‏ چین‌خوردگی و جابه‌جایی رسوبات نشوژن گاهی 
جوان‌تر است. در اثر فعالیت این گسل واحدهایی از کنگلومرا؛ مارن گچدار ماسه سنگ خاکستری» مارن‌های با Ol‏ 
Sal al‏ و همچنین مارن خاکستری گچدار و نمک‌دار در کنار هم قرار گرفته‌اند. با بررسی معیار فاصله از گسل به 


این نتیجه دست‌يافتيم که طبق نقشه پهنه‌بندی خطر فرونشست به دست آمده از روش مارکوس. پهنه‌های با احتمال 


سال یازدهم پهنه بندی گسترهٌ خطر فرونشست زمین در دشت سراب .... yw‏ 


خطر زیاد و بسیار زیاد در فاصله صفر تا ۸۰۰۰ متری» از گسل واقع شده‌اند و بر اساس نتایج حاصل از روش 
کوداس. مناطق پرخطر و بسیار پرحطر در فاصله صفر تا ۱۲۰۰ متری» از گسل قرار دارند. 

عامل نزدیکی به رودخانه به‌عنوان عامل مزاحم در ایجاد فرونشست مطرح می‌باشد و با فاصله گرفتن از رود امتیاز 
بیشتری به زمین‌های منطقه تعلق می‌گیرد. هر چه از آبراهه‌ها فاصله می‌گیریم نیاز استفاده از آب‌های زیرزمینی به دلیل 
عدم وجود و استفاده از آب‌های سطحی بالا می‌رود و در پی استفاده نادرست از آب‌های زیرزمینی و افت سطح 
آب‌های زیرزمینی» میزان احتمالل فرونشست افزایش می‌یابد. رودخانه‌های اصلی جاری در دشت سراب؛ بیوک‌چای, 
تاجیارچای» رازلیق چای» وانق‌چای می‌باشد و منابع تغذیه دشت عمدتاً آب‌های نفوذیافته از جریان‌ه ای سطحی 
ارتفاعات حاشیه دشت و همچنین نفوذ آب حاصل از باران بر سطح دشت است. هر چه از آبراهه‌ها فاصله می‌گیریم 
نیاز استفاده از آب‌های زیرزمینی به دلیل عدم وجود و استفاده از آب‌های سطحی Ub‏ می‌رود و در پی استفاده نادرست 
از آب‌های زیرزمینی و افت سطح آب‌های زیرزمینی» ميزان احتمال فرونشست افزایش می‌یابد. به‌طور AS‏ فاصله از 
آپراهه‌ها در محدوده مطالعاتی بین صفر تا 2۰۰۰ متر می‌باشد و با انطباق نقشه پهنه‌بندی خطر فرونشست حاصل از 
هر دو روش مورد استفاده در مطالعه حاض با نقشه معیار فاصله از رودخانه؛ می‌توان گفت که به‌طور عمده مناطق 


دارای احتمال be‏ در فاصله صفر تا ۰ متری از رودخانه قرار دارند. 


-٤‏ جمع‌بندی 
در این پژوهش با استفاده از الگوریتم‌های تصمیم‌گیری چند معیاره CODAS‏ و MARCOS‏ به ارزیابی و 
پهنه‌بندی as‏ فرونشست دشت سراب. پرداخته شده است. بررسی‌های انجام‌گرفته در این تحلیل نشان داد که به 
ترتیب؛ افت سطح eol‏ کاربری اراضی؛ شیب و زمین‌شناسی در وقوع فرونشست محدوده مطالعاتی؛ بیشترین اهمیت 
را دادند و می‌توان گفت که مهم‌ترین عامل اصلی دخیل در افزایش مقدار و پتانسیل فروند ت دشت سراب» 


بهره‌برداری بی‌رویه از آب زیرزمینی می‌باشد. در راستای تأیید نتیجه حاصله بایستی اشاره داشت که شرکت آب 


منطقه‌ای آذربایجان شرقی نیز اعلام os S‏ است که سالانه ۱۵ میلیون مترمکعب آب از چاه‌های مجاز سراب اضافه 
برداشت می‌شود و همچنین» ناموری و صدریکیا (۱۳۹۹)» صدری‌کیا OE)‏ در پژوهشی که در رابطه با بررسی 
فرونشست دشت سراب انجام داده‌انده به این نتیجه دست adh‏ بودند که عدم کنترل و مدیریت مناسب منابع آب‌های 
زیرزمینی و افت سطح آب زیرزمینی به افزایش میزان فرونشست در دشت سراب منجر شده است. لذا نتیجۀ پژوهش 
حاضر با نتایج این مطالعات cd‏ هم‌خوانی دارد. 

به‌علاو» خروجی حاصل از هر دو روش CODAS MARCOS‏ نشان‌دهندة توان بالای دشت سراب از 
لحاظ خطر فرونشست می‌باشد و با توجه به این امر که فرونشست زمین می‌تواند خسارت‌های جبران ناپذیر مالی و 


جانی به دنبال داشته باشد و به بسیاری از سازه‌های سطحی و زیرسطحی در مناطق شهری و Ao io.‏ آن» در محدوده 


\w‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


موردمطالعه سیب tly‏ لازم است تا سازمان‌های متولی فرونشست زمین» مثل؛ سازمان زمین‌شناسی و اکتشافات 
معدنی. شرکت آب منطقه‌ای, جهاد کشاورزی, حفاظت از محیط‌زیست» منابع طبیعی و امور آبخیزداری و 
شهرداری‌هاء با انجام اقدامات مسئولانه و هم‌افزایی لازم مانع کاهش ذخایر آبی و فرونشست‌های خطرآفرین برای 
حیات در منطقه انجام دهند. ازاین‌رو بایستی راهکارهای بلندمدت و کوتاه مدت. به‌صورت اصلاح روش‌های مدیریت 


منابع آب. جلوگیری از برداشت سفره‌های آب زیرزمینی» توقیف تجهیزات حفاری غیرمجاز, ساماندهی شرکت‌های 


حفاری» پایش مستمر مناطق مستعد فروند ت» ایجاد سدهای زیرزمینی به‌منظور افزایش تراز آب و کاهش 
"Tc‏ ت» مهار آب‌های سطح و اجرای طرح‌های تزریق مصنوعی آب به سفره‌های زیرزمینی منطقه (طبق گفتة 


رئیس اداره منابع آب شهرستان سراب که با توجه به آورده رودخانه‌های محلی» می‌توان سالانه حدود ۳ میلیون 
مترمکعب آب به این سفره‌ها تزریق OS‏ تهیه و نصسب کنتور هوش‌مند و کنترل برداشست ol‏ آموزش عمومی 
کشاورزان محلی. کنترل شدید مصرف آب و تغیبر الگوی مصرف. مورد توجه مسئولان ذی‌ربط و مربوطه و همچنین 
ساکنان محلی قرار گیرد. 

درنهایت قابل ذکر است. همان‌گونه که می‌دانیم» هر یک از روش‌های dallas‏ مخاطرات محیطی ازجمله 


برخوردار است به شرایط طبیعی و انسانی محدوده بستگی دارد. انتظار می‌رود. نتایج حاصل از پژوهش حاضر با 
معرفی روش‌های تحلیل چند معیاره MARCOS‏ و ۵6 که تاکنون در هیچ مطالعه‌ای با هدف shig‏ 
فرونشست مورد استفاده قرار نگرفته بوده‌انده در امر انتخاب روش مناسب بررسی و پهنه‌بندی خطر فرونشست مورد 


توجه پژوهشگران قرار گیرد. 
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